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ABSTRACT 

We study the thermal structure of neutron stars with magnetized envelopes composed of accreted 
material, using updated thermal conductivities of plasmas in quantizing magnetic fields, as well as 
equation of state and radiative opacities for partially ionized hydrogen in strong magnetic fields. The 
relation between the internal and local surface temperatures is calculated and fitted by an analytic 
function of the internal temperature, magnetic field strength, angle between the field lines and the 
normal to the surface, surface gravity, and the mass of the accreted material. The luminosity of a 
neutron star with a dipole magnetic field is calculated for various values of the accreted mass, internal 
temperature, and magnetic field strength. Using these results, we simulate cooling of superfluid 
neutron stars with magnetized accreted envelopes. We consider slow and fast cooling regimes, paying 
special attention to very slow cooling of low-mass superfluid neutron stars. In the latter case, the 
cooling is strongly affected by the combined effect of magnetized accreted envelopes and neutron 
superfluidity in the stellar crust. Our results are important for interpretation of observations of 
isolated neutron stars hottest for their age, such as RX J0822— 43 and PSR B1055— 52. 
Subject headings: dense matter — magnetic fields — stars: individual (RX J0822. 0—4300, PSR 
B1055— 52) — stars: neutron 



1. INTRODUCTION 

Cooling of neutron stars (NSs) depends on the prop- 
erties of dense matter in their crusts and cores. The- 
oretical predictions of these properties depend on a 
model of strong interaction and a many-body theory 
employed. Therefore, one can test microscopic mod- 
els of dense matter by comparing the results of sim- 
ulations of NS cooling; with observations of the stella r 
thermal emission (e.g., IPagellT997l Il99l IPethicldli99l . 
The history of NS cooling theory is re viewed, e.g., by 
lYakovlev. Levenfish. fc Shibanovl l)1999l) ; the main cool- 
ing reg ulators are discussed, e.g., by lYakovlev et al.l 
(2002bJ. Practical implementation of this method is re- 
stricted by the accuracy of observational data and by the 
quality of the physics input used in calculations. 

It is conventional (e.g., 

iGudmundsson. Pethick. &: Epsteinlll983|) to separate the 
outer heat-blanketing envelope from the stellar interior. 
The blanketing envelope can be treated separately 
in the plane-parallel qua si-Newtonian approximation 
( Gud mundsson et al.lll983f) , which simplifies the cooling 
calculations. 

The NS cooling is strongly affected by the relation be- 
tween the heat flux density F and the temperature at 
the inner boundary of the blanketing envelope, Ti>. The 
flux density is constant throughout a given local part of 
the envelope. It is conventionally written as F = crT 4 , 
where T s is the local effective surface temperature, and a 
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is the Stefan-Boltzmann constant. The T\,—T s relation is 
determined by the equation of state (EOS) and thermal 
conductivity of matter in the heat-blanketing envelope. 

For nonmagnetized stellar envelopes com- 
posed of iron, thi s relationship wa s thor - 
oughly studied b v IGudmundsson et alJ \ 19831). 
Cha brier. Potekhin. k Yakovlev ! il997l) and 

IPotekhin. Chabrier. fc Yakovlevl l|1997fl (hereafter 
Paper I), reconsidered the problem and extended the 
results in two aspects. First, advanced theoretical 
data on EOS and thermal conductivity allowed the 
authors to study colder NSs, with T s down to 5 x 10 4 K. 
Second, the authors considered the blanketing envelopes 
composed not only of iron but also of light (accreted) 
elements. They found that a small amount of accreted 
matter (with mass 10~ 16 M Q < AM < 10" 7 M Q ) result 
in higher T s at a given Tt,. As a consequence, the stellar 
luminosity is higher at the early (neutrino-dominated) 
cooling stage (at stellar ages t < 10 4 — 10 5 yr) and lower 
at the subsequent photon stage. 

Recently the ther mal structure of accreted en v elopes 
has been studied by iBrown. Bildsten. fc Changl (2002) 
taking into account t hat hydrogen burn ing in accreted 
matter may proceed llSchatz et al.1 1200 i]) far beyond Fe 
(up to Te) via the rapid prot on capture process. 

It was realized long ago ijTsuruta et alJ [1972) that 
strong magnetic fields can greatly affect the surface tem- 
perature and cooling of NSs. The thermal structure 
of NS envelopes with radial magnetic fields ( normal to 
the surface) was studied by iVan Rineil l)1988f ) (also see 
I Van Riper! 119881 for references to the earlier work). His 
principal conclusion was that the field reduces the ther- 
mal insulation of the blanketing envelope due to Landau 
quantization of electron motion. 

The thermal structure of the envelope with magnetic 
fie lds normal and ta ngen tial to the surf ace was analyzed 
bv lHernauistl ij 19 85ft and lSchaaj l|1990aft . The tangential 
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field increases the thermal insulation of the envelope, be- 
cause the Larmor rotation of the electrons reduces the 
transverse electron thermal conductivity. 

The case of arbitrary angle 9 between the field 
lines and the normal to the surface was studied by 
IGreenstein fc Hartkei <| 198.11) in the approximation of con- 
stant (density and temperature independent) longitudi- 
nal and transverse thermal conductivities. The authors 
proposed a simple formula which expresses T s at arbi- 
trary 9 through two values of T s calculated at 9 = and 
90° . Th e ca se of arbitrary 9 was studied a lso bv iSchaal 
(199013) and IHevl fc Hernmiisti GMIIlLl). 

iPotekhin fc Yakovlev! (|200lD (hereafter Paper II) re- 
considered the thermal structure of the blanketing en- 
velo pes for any 9, u sing improved thermal conductivi- 
ties IPotekhinl 1999). In agreement with an earlier con- 
jecture of iHernauistl 1119851) and simplif i ed tre atments of 
IPaeel (|1995() and lShibanov fe Yakovlev! |l99^ . Paper II 
demonstrated that the dipole magnetic field (unlike the 
radial one) does not necessarily increase the total stellar 
luminosity L 7 at a given T^. On the contrary, the field 
B ~ 10 n -10 13 G lowers L 7 , and only the fields B > 10 14 
G significantly increase it. 

Early simulations of cooling of magnetized NSs were 
performed assuming the radial magnetic field everywhere 
over the stellar surf ace fe.g.. iNomoto fc Tsurutalll987t 
iTsuruta et all Il972t iVan Rioerl 119911) . Since the ra- 
dial magnetic field reduces the thermal insulation, these 
models predicted the increase of L 7 by the mag netic 
field at the neutrino coolin g stage. iPagd l) 19951) and 
iShibanov fc Yakovlev! (1996) simulated cooling of NSs 
with dipole magnetic fields, using the Greenstein-Hartke 
formula and the T^-T s relations at the ma gnetic pole 
and equator taken from the p reviou s work ( Hcrnauist 
T^IS^haajll990atlVan Ripeiill988|) . iHevl fc Hernouistl 



1997a b, 1998, 2001) proposed simplified models of cool- 
ing magnetized NSs including the cases of ultrahigh sur- 
face magnetic fields, B ~ 10 15 -10 16 G. Recently, the im- 
proved Tb-T s relation derived in Paper II has been used 
for simulating th e cooling of NSs with dipole magnetic 
fields (Paper I I; Kaminker. Yakov lev. fc Gnedinl 120021: 
lYakovlev et al.H2002bl) . 

Most of the above-cited work on magnetized NSs has 
been focused on NSs wit h iron envelopes, except that 
IHevl fc Hernauisti l)1997a|) have used simplified analytic 
models to calculate cooling of ultramagnetized (B ~ 
10 15 -10 16 G) NSs which may possess accreted envelopes. 

In this paper, we consider accreted NS envelopes 
(as in Paper I for B = 0), but take into account 
strong magnetic fields B directed at various angles 
9 (as in Paper II for iron envelopes). We use the 
modern electron conductivities of magnetized envelopes 
(jPotekhinl I1999J) . taking into account the effects of fi- 
nite size of atomic nuclei, appropriate for the inner en- 
velopes (jGnedin. Yakovlev. fc Potekhinl 120011) . as well 
as the radiative opacities and EOS of strongly mag- 
netized, partially ionized h ydrogen in NS atmosphere 
l)Potekhin fc Chabrierl l2003) . We derive an analytic ap- 
proximation for the surface temperature. We calculate 
T e g for a dipole magnetic field and simulate NS cooling 
to demonstrate the effects of magnetic fields and accreted 
envelopes on thermal evolution of NSs. 
In f|21we give an overview of the main effects of strong 



magnetic fields on heat conduction in a plasma, and in- 
troduce basic definitions. In fj^jwe describe the physics 
input used in our simulations. The thermal structure of 
NS envelopes with magnetic fields is discussed in <2] NS 
cooling simulations are presented in ^ The conclusions 
are formulated in ^ 

2. THE EFFECTS OF STRONG MAGNETIC FIELDS ON HEAT 
CONDUCTION IN NEUTRON STAR ENVELOPES 

We call a magnetic field strong if the electron cyclotron 
energy fku c = heB/m e c exceeds 1 a.u. — i.e., the field 
strength B is higher than mice 3 /h 3 = 2.3505 x 10 9 G, 
where m G is the electron mass, e the elementary charge, 
and c the speed of light. We call the field super strong if 
Hlu c > m c c 2 , that is B > m 2 c c 3 /eh = 4.414 x 10 13 G. 

The field is called strongly quantizing if it sets almost 
all electrons on the ground Landau level. The latter case 
takes place at temperature T <C Tg and density p < ps, 
where 

PB=m u n B (A)/(Z) « 7045B 3 2 /2 ((A)/(Z)) g cm- 3 {l) 

Tb = hujc/k-Q w 1.343 x 10 8 B\2 K. (2) 

Here, m u = 1.66054 x 10 -24 g is the atomic mass unit, 
(A) and (Z) are the mean ion mass and charge num- 
bers, B 12 = B/(10 12 G), n B = l/(7r 2 V2a 3 n ) is the elec- 
tron number density at which the electron Fermi en- 
ergy reaches the first excited Landau level, and a m = 
(Hc/eB) 1 / 2 is the quantum magnetic length. 

As a rule, at p > ps and not too high T, the elec- 
tron thermodynamic and kinetic functions oscillate, with 
increasing p, around their values calculated neglecting 
quantization of electron motion in Landau levels. If 
T > Tb or p ^> ps, the field can be treated as classi- 
cal (nonquantizing). 

The increase of the atomic binding energies in the 
strong magnetic field tends to lower the ionization de- 
Tee, a s first guessed bv lCohen. Lodenauai. fc Rudermanl 
1970). Hence, there can be a significant amount of 
bound species in a highly magnetized atmosphere, even 
if it were negligibly small at the same temperature in the 
zero-field case. This idea has recen tly been confirmed by 
detailed N S atmosphere models (IPotekh in fc Ohabrieil 
2001 IPotekhin. Chabrier. fc Shibanovl Il999ai 
Raiagppal. Romani. fc Milleilll997ft . 



The effects of magnetic fields on the kinetic proper- 
ties of the outer NS layers ha ve be en reviewed, e.g., by 
Yakovlev fc Kaminkerl l)1994|) and iVentura fc Pot ekhin 
( 2001). A strongly quantizing magnetic field reduces 
the m ean radiative opacities (e.g.. iSilant'ev fc Yakovlev! 
1980); therefore the b ottom of the photosphere shifts to 
highe r densities (e.g.. lLai fc Saloeterlll997HPavlov et alJ 
1995). 

A nonquantizing magnetic field does not affect thermo- 
dynamic functions of the plasma. However, it does affect 
the electron heat conduction, unless the magnetization 
parameter 
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(3) 



is small. Here, tu g = ui c / + x 2 is the electron gyrofre- 
quency, 

x x = ft(37r 2 n ) 1/3 /(m c c) w 1.0088 {p 6 Z/A) 1/3 (4) 
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is the relativity parameter, p$ = p/(10 6 g cm -3 ), and r 
is an effective relaxation time. In a degenerate Coulomb 
plasma with a nonquantizing magnetic field, 



3ttH 3 



5.700 x 10" 



4Zm e e 4 A ^1 + x 2 ZA y/T 



(5) 



where A is an effective Coulomb logarithm. One has 
A ~ 1 in the ion liquid, and A is roug hly proportional to 
T in the classical Coulomb solid ("e.g , ..lYakoy ley_fc U rnir " 
1980). However, according to iBaiko et al . (1998), the 
correlation effects smooth the dependence of A on T or 
p near the melting point and introduce deviations from 
this proportionality. 

In a magnetic field, one should generally introduce two 
different relaxation times (tii and tj_), determined by 
two Coulomb logarithms (A|| and Aj_), corresponding to 
electron transport parallel and perpendicular to the field 
lines. Analytic fits for A| | and Aj_ at any B have been 
constructed bv lPotekhinl |l299). 

As seen from Eqs. (0 and ©, the magnetization pa- 
rameter is large in the outer NS envelope at B > 10 11 G, 
typical for many isolated NSs. Moreover, according to 
Eqs. and @ the magnetic field can be strongly quan- 
tizing in the outermost part of the envelope. Therefore, 
the field can greatly affect the heat conduction and the 
thermal structure of the NS envelope. In this case, the 
surface temperature T s may be nonuniform, depending 
on the magnetic field geometry. Then, we introduce the 
effective temperature of the star T c g defined by 



47rcri? 2 T c 4 ff = L 7 = / FdE = a I T 4 d£ 



(6) 



where R is the stellar (circumferential) radius, and dS 
is the surface element. The quantities T s , T c ff, and L 7 
refer to a local reference frame at the NS surface. The 
redshifted ("apparent") quantities as detected by a dis- 
tant observer are l)Thornelll9"77|) : T* 

(1 



1 cS 



T s y/1 - r g /R, 

T c tf a/1 - Tg/R, and = (1 - r g /R) i 7 , where 
r g = 2GM/c 2 = 2.95(M/M Q ) km is the Schwarzschild 
radius defined by the total gravitational NS mass M, and 
G is the gravitational constant. 

3. PHYSICS INPUT 

Much of the physics input used in our work is the same 
as in Papers I and II. We outline these cases and refer 
the reader to Papers I and II for details. 

3.1. Basic Equations 

We choose the boundary between the internal NS re- 
gion and the heat-blanketing envelope at the neutron- 
drip density pb = 4 x 10 11 g cm -3 (at a radius Rb a 
few hundred meters und er the surface). Although p b = 
10 10 g cm~ 3 adopted bv iGudmundsson et alJ l)1983[) en- 
sures the necessary condition R—Rb <C R with higher ac- 
curacy, the present choice better conforms to the require- 
ment that T does not vary over the boundary, if a strong 
magnetic field is present (Paper II). In addition, we show 
that this choice provides a better accuracy to the T s -Tb 
relation in the presence of an accreted envelope. On the 
other hand, the increase of pb increases the thermal re- 
laxation time of the blanketing envelope, thus hamper- 
ing a study of rapid variability of the thermal emission 



(over time scales < 10 yr fo r pb ~ 4 x 10 11 g cm~ 3 , 
HJshomirskv fc RutledgelEOQlfl with a cooling code. 

The thermal structure of the blanketing envelope is 
studied in the stationary, local plane-parallel approxi- 
mation, assuming that a scale of temperature variation 
over the surface is much larger than the thickness of the 
blanketing envelope. This leads to the one- dimensional 
approximation for the heat diffusion equation: 



dT 

dz' 



16erT 3 
3Kp ' 



(7) 



where k is an effective thermal conductivity along the 
normal to the surface, K is the mean opacity, and z = 
(R — r) (1 — rg/R)^ 1 / 2 <C -R is the local proper depth in 
the envelope. 

Equation @ can be reduced fe.g-. lVan Riperlll988f) to 



dlogT 
dlogP 



3 PK T 4 
16 ~g~ T*'' 



(8) 



where P is the pressure. We integrate Eq. (JHJ inwards 
from the radiative surface, where T = T s , using the 
same algorithm as in Paper I. We place the radiative 
surface at the Rosseland optical depth equal to 2/3. 5 
Another boundary condition would imply another inte- 
gration constant for the integral of Eq. ©, which be- 
comes negligi ble in the inner part of the envelope, where 
T > T s re.g.. lVentura fc Potekhinll2001fl . Therefore the 
choice of the boundary condition does not noticeably af- 
fect the heat flux and the T s -Tb relation. 

The validity of the one-dimensional approximation can 
be checked by two-dimensional simulation of the heat 
transport in the bl anketing envelo pe. Such a simulation 
was attempted by ISchaaJ |l990bj) for a homogeneously 
magnetized NS under many simplified assumptions. The 
heat conduction from hotter to cooler zones along the 
surface or possible meridional and convective motions 
can smooth the temperature variations over the NS sur- 
face. Nevertheless, the one-dimensional approximation 
seems to be sufficient to simulate NS cooling. 

The thermal conductivity tensor of a magnetized 
plasma is anisotropic. It is characterized by the con- 
ductivities parallel (/tii) and perpendicular (k±) to the 
field, and by the off-diagonal (Hall) component. In the 
plane-parallel approximation, Eq. (JJJ contains the effec- 
tive thermal conductivity 



K = Ku cos 



k± sin 



(9) 



Assuming the di pole field, we use the general- 
relativistic formulas (Gin zburg fc Ozernovlll96 4) 



B (x) = Bp \J cos 2 x + a 2 sin 2 x, 
tan0 = atanx, 
_ (1 - x) ln(l - x) + x - 0.5x 2 
[ln(l - x) + x + 0.5x 2 ]VT~^ x' 



(10) 
(11) 
(12) 



where B p is the field strength at the magnetic pole, x 1S 
the polar angle, and x = r g / R. 

5 Assuming the local thermodynamic equilibrium, this depth 
follows from the Milne— Eddington formula, which is fairly accurate 
as shown by Koureanoff (1948). 
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Fig. 1. — Temperature profiles in a nonaccreted, nonmagnetized 
envelope of a "canonical" NS (Af = IAMq, R = 10 km) for several 
values of T B (marked by logT s ). Solid lines: calculation with the 
improved electron conductivities; dashed lines: Paper I. The dotted 
line shows the liquid/solid boundary. Inset shows a temperature 
profile near this boundary in more detail. 




logp (g cm" 3 ) 

Fig. 2. — Temperature profiles in a nonmagnetized envelope of 
the canonical NS containing maximum accreted mass (solid lines), 
compared with the iron envelope (dot-dashed lines). The dotted 
lines correspond to interfaces between the layers of different chem- 
ical elements. The inset demonstrates the temperature increase in 
the hot envelope behind the deepest accreted layer. 



3.2. Equation of state 

A nonaccreted heat-blanketing envelope is assumed 
to be composed of iron, which can be partially ion- 
ized at p < 10 6 g cm~ 3 . Ions in the envelope form 
either Coulomb liquid or Coulomb crystal. The melt- 
ing of the crystal occurs at T = T m w 175 (e.g., 
iPotekhin fc Chabriedl2T)0TL where T = (Ze) 2 /k B T ai is 
the ion coupling parameter, a\ — {Aim- l /'i)^ 1 ^ is the ion 
sphere radius, and n\ is the ion number density. The 
melting curve T = T m is shown in Fig. ^ by the dot- 
ted line. The figure demonstrates the thermal struc- 
ture of a nonmagnetic iron envelope. The value of the 
surface gravity g — 2.43 x 10 14 cm s -2 chosen in this 
example corresponds to the "canonical" NS model with 
M = 1.4 Mq and R = 10 km. At B = 0, as in Paper I, w e 
use the OPAL EOS l|Rogers. Swenson. fc Iglesiasl [l996). 
extrapolated beyond the available tables using an effec- 
tive charge number Z e g. In the strong magnetic field 
w e employ the finite-temp erature Thomas-Fermi EOS 
of Thorolfsson et al. (1998), 6 with Z c g evaluated as in 
Paper II. 

As in Paper I, the accreted envelope is modeled by a 
sequence of layers of H, He, C, O, and Fe; and Fe is con- 
sidered as the end point of nuclear transformations in the 
heat-blanketing envelopes. An example of the thermal 
structure of a nonmagnetic accreted envelope is shown 
in Fig. |3 The boundaries between the layers are de- 
termined by the conditions of thermo- and pycnonuclear 

6 Available at http://www.raunvis.hi.is/~ath/TFBT/ 



burning of the selected elements and by AM, the mass of 
light elements (from H to O). These approximate bound- 
aries are shown in Fig. [21 by dotted lines. In that exam- 
ple, the envelope is fully accreted: light elements reach 
the density ~ 10 10 g cm~ 3 , where pycnonuclear burning 
of oxygen becomes efficient. Compared to Paper I, we 
have improved the position of the C/O boundary , takin g 
into account the results of iSahrling fc Chabrierl fl998). 
Their carbon ignition curve is approximately reproduced 
by the formula 



T 



5.2 x 10 s K 



{l + [0.2 1n(p pyc /p)]-"-7} 



0.710.07 ' 



(13) 



at p < p pyc w 5.5 x 10 9 g cm 3 . The latter improvement, 
however, has almost no effect on 



As demonstrated recently bv iSchatz et alJ 1)200 1|) . an 
explosive or steady-state burning of hydrogen in the sur- 
face layers of accreting NSs may be strongly affected by 
the rapid proton (rp) capture process and extend to the 
elements like Sn, Sb, and Te, much heavier than Fe. We 
do not consider this possibility here but intend to analyze 
it in a separate publication. 

The outermost accreted layer is assumed to consist of 
hydrogen. Even a hydrogen atmosphere of a "warm" 
NS can be partially ionized in a strong magnetic field. 
Therefore we use the EOS of partially ionized hydro- 
gen in strong field s, derived by IPotek hin et alJ £l999a) 
and tabulated bv IPotekhin fc Chabrier (2003) at 11.9 < 
log 10 B < 13.5 (where B is in G). Beyond the tabu- 
lated range we employ the model of a fully ionized ideal 
electron- ion plasma. As demonstrated below (Fig. |3J), 
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Fig. 3. — Temperature profiles (upper panels) in the accreted envelope (AM/M = 1.5 X 10 7 ) of the canonical NS with B = 
panels) and B = 10 13 ' 5 G (right), for T B = 2 X 10 5 K (the curves which start in the atmosphere at lower T) and T s = 10 6 K. 
the lower curve corresponds to the radial magnetic fi eld (8 = 0, ||) and the upper to the tangential one (9 = 90°, _L). The solid lines are 
obtained using the hydrogen EOS and the opacities of Potckhin & Chabricr (2003) for the partially ionized atmosphere; the corresponding 
mean ion charge number is plotted on the lower panels. The dashed profiles are calculated using the EOS of fully ionized ideal gas. The 
two arrows marked 'ps' indicate the density below which the field is strongly quantizing for hydrogen (left arrow) or helium (right arrow 
on each panel). The arrows marked '1' and '2' indicate filling of the first and the second electron Landau levels, respectively (ps and 
(2 + v2) PSi calculated for H on the left and for 12 C on the right panel). 



the inaccuracy introduced by the latter approximation 
can be significant for cold NSs (7b < 10 7 K) with super- 
strong magnetic fields (B > 10 13 5 G). The fully ionized 
ideal electron-ion plasma model is used also for He, C, 
and O, which occupy deeper layers. 

3.3. Opacities 

The heat is carried through the NS envelope mainly 
by electrons at relatively high densities and by photons 
near the surface. Hence, 

K = Kl + Ke , K- 1 = K~ 1 +K- 1 , (14) 

where K r , n c and K T , K c denote the radiative (r) and 
electron (e) components of the thermal conductivity k 
and opacity K. 

Typically, the radiative conduction dominates (ft r > 
Ac e ) in the outermost nondegenerate NS layers, whereas 
the electron conduction dominates (k c > n r ) in deeper, 
moderately or strongly degenerate layers. We will see, 
however, that in the superstrong magnetic fields the ra- 
diative conduction can be important up to higher den- 
sities. The Tb-T s relation mostly depends on the con- 
ductivities in the sensitiv ity strip on the p—T plane 
ijGudmundsson et al.lll983|) . where k c ~ « r . The T(p) 
profiles flatten in the inner NS zone beyond this strip 
(e.g., see Figs. Ql EJ. 

Following Paper II, we us e the updated e l ectron con- 
ductivities re e , presented bv iPotekhin et all (|1999bf) for 



B = and lPotekhinl (|1999^ for B ^ 0. 7 Unlike the pre- 
vious expressions for K e , the updated results take into 
account multiphonon absorption and emission processes 
in Coulomb crystals and incipient long-range order in 
strongly cou pled Coulomb liquids of ions. As argued by 
iBaiko et all l|1998j) . these processes are important near 
the melting, 0.3 T m < V < 3T m , and they almost remove 
the jump of K e at the liquid-solid interface. 

The effect of the conductivity update on the ther- 
mal structure at B — is illustrated by Fig. 2] It 
is marginally significant at low Tb and insignificant at 
T b > 10 s K. 

The inset in Fig. [21 zooms the innermost parts of the 
upper (hottest) profiles for the nonaccreted and accreted 
envelopes. For the accreted envelope, the profile notice- 
ably mounts at p > 10 10 g cm~ 3 , because the crossing 
of the O /Fe interface is accompanied by a decrease of k c 
due to the jump of Z . The decrease of K e has smaller 
impact on the temperature profiles at lower T. 

Radiative opacities of hydrogen at moderately strong 
magne tic field are taken from IPotekhin fe Chabrierl 
(2003). These opacities are contributed by free- free, 
bound-free, and bound-bound transitions, and electron 
scattering. They are calculated taking into account de- 
tailed ionization balance in strong magnetic fields. Fig- 
ure shows the thermal structure of accreted NS en- 

7 Available at http://www.ioffe.ru/astro/conduct/ 
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velopes at two values of T s and two values of B, for 
the magnetic fields radial and tangential to the NS sur- 
face. S olid lines are obtained using the EOS and opac- 
ities of iPotekhin fc Chabrierl l|2003D . while dashed lines 
correspond to the ideal fully ionized plasma model. At 
T s = 10 6 K the latter model appears to be satisfactory, 
whereas at lower T s there are appreciable differences. 

The difference at the lowest densities (in the atmo- 
sphere) is produced by the contribution of bound species 
(atoms) in the EOS and opacities. The lower panels of 
Fig. OH show the mean ion charge (Z) along the lower 
(T s = 2 x 10 5 K) profiles, confirming that (Z) is consid- 
erably smaller than 1 in the atmosphere. At this rela- 
tively low temperature, the ionization proceeds smoothly 
for B = 10 12 G, but rather sharply (via pressure ioniza- 
tion) for B = lO 13 5 G. Simultaneously with the pressure 
ionization, the temperature steeply climbs along the pro- 
files on the upper right panel, because a sharp increase 
of P (due to the increasing number of free particles) con- 
tributes to the right-hand side of Eq. (JHJ . 

The difference between the ideal-gas (dashed) and ac- 
curate (solid) profiles at higher (subphotospheric) densi- 
ties arises from the Coulomb interaction. This interac- 
tion gives a negative contribution to P, thus decreasing 
the right-hand side of Eq. (JSJ) . 

We adopt the same radiative opacities of iron as in Pa- 
pers I and 11 for zero and strong B, r espectively. In the 
first c ase, we use the OPAL opacities (Iglesia s fc Rogers! 
1996). For strong B, we use the free-free and scattering 
opacity fits from Pap er II. However, in the presen t paper 
(as in Paper I and in IPotekhin fc Chabrierll200^ we as- 
sume that radiation does not propagate at frequencies be- 
low the electron plasma frequency w p i = (4ire 2 n c /m c ) 1 / 2 ; 
thus we cut off the integration over photon frequencies u> 
in the Rosseland mean opacities at uj < ui p \. We found 



that, under this assumption, the fit expressions for K T 
given in Paper II should be multiplied by a correction 
factor w exp{0.005[ln(l + 1.5huj pl /k B T)} 6 }, which effec- 
tively eliminates the radiative transport at large densi- 
ties, where hu) v \ w 28.8 keV^p 6 Z/A > k B T. 

For He, C, and O, we also employ the opacity fits from 
Paper II with the above plasma-frequency correction. 
The effect of this correction on the thermal structure of 
the accreted envelope is illustrated in Fig. 0] It is unim- 
portant at moderately strong magnetic fields, but quite 
significant at superstrong fields, because of the high pho- 
tosphere densities in the latter case. The effect becomes 
less pronounced at higher T s . Actual ly some plasma 
waves can propagate at u> < lu p i (e.g., iGinzbura 11970(1 
and carry the heat; thus we expect that realistic tem- 
perature profiles should lie between the solid and dashed 
curves in Fig. 0] 

4. THERMAL STRUCTURE 

4.1. Temperature profiles 

We integrated Eq. JHJ) for various magnetic field 
strengths B, inclination angles 9, surface temperatures 
T s , and accreted masses AM, using the envelope mod- 
els described in the previous section. We performed 
the calculations at two values of the surface gravity, 
g = 10 14 cm s -2 and 2.43 x 10 14 cm s~ 2 , and checked 
th at the approximate scaling r elation T s cx g 1 ^ 4 obtained 
bv lGudmundsson et al.1 1)1983|) and confirmed in Papers I 
and II, holds also in the present case. 

As clear from the discussion in SJ31 °ur EOS and opacity 
models may be crude at B > 10 14 G. Further improve- 
ments of physics input are required in superstrong fields. 

Some examples of calculated temperature profiles at 
several values of T s are shown in Figs. and discussed 
above. Figure [S] displays the profiles at two values of B 
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Fig. 5. — Temperature profiles in the accreted envelope for the 
same stellar parameters as in Fig. but at B = 10 12 G (left 
panels) and B = 10 14 G (right), for T b = 10 8 K (lower panels) and 
10 9 K (upper ones), at 8 = 90° (steepest lines), 84° (middle), and 
(the least steep lines). The dot-dashed lines correspond to iron 
envelopes at 8 = 84° . 



and three values of 9, for AM = 2x10 7 M Q (solid lines) 
and (dot-dashed lines), and for two values of Tb (10 
K and 10 9 K). 

4.2. Relation between internal and effective 
temperatures 

Using the above results we have calculated the Tt,-T s 
relation for a number of input parameters. We confirm 
the conclusion of Paper I that the relation is rather insen- 
sitive to the details of the distribution of different chem- 
ical elements within the accreted envelope but depends 
mainly on AM, the total mass of the elements from H 
to O. In particular, all hydrogen can be replaced by He 
leaving the Tb~T s relation almost unchanged. 

We h ave compared our Th —T s relation with that calcu- 
lated bv lBrown et alJ <)2002j) (their Fig. 4) for nonmagne- 
tized accreted envelope composed of H, He, and Fe. The 
agreement is quite satisfactory. 

We have fitted our numerical Tb-T s relation calculated 
for magnetized accreted envelopes by analytic formulas 
(presented in the Appendix). In the limits of B = 
and/or AM — these formulas do not exactly repro- 
duce the fits obtained in Papers I and II. The differences 
reflect the improvements in the physics input and enve- 
lope models, discussed in For instance, in calcula- 
tions shown in Fig. the boundary condition for the 
integration inside the star was determined using the fit- 
ting formulas for T s , presented in the Appendix. The 
good convergence of the profiles toward the desired Tb 
confirms the accuracy of our fits. 



Figure El illustrates the Tb~T s relation for the fully ac- 
creted envelope of the canonical NS without magnetic 
field, with moderately strong field B = 10 12 G, and with 
superstrong field B = 10 15 G, for two field geometries 
- normal (9 = 0) and tangential (9 = 90°) to the sur- 
face. The lines show the fit, and the symbols (dots and 
triangles) show the numerical results. The nonmagnetic 
fit of Paper I is plotted by the dashed line. Its deviation 
from the present results (solid line) at high temperatures 
is mainly explained by the change of ph ( jSH) - 

Figurc[7|illustrates the Tb-T s relation for the canonical 
NS with fully, partly, and nonaccreted envelopes at 9 = 
and 90°. The longitudinal heat transport (9 = 0) is en- 
hanced by the accreted envelope and/or by the strong 
magnetic field, thus increasing T s . However, the trans- 
verse transport (lower lines), which is reduced by the 
strong magnetic field, is additionally reduced by the ac- 
creted envelope, further decreasing T s . Actually, in lon- 
gitudinal and transverse cases the lowering of Z increases 
the effective electron relaxation time r. However, at 
WgT 3> 1 the transverse electron conductivity (contrary 
to the longitudinal conductivity or the condu ctivity at 
B = 0) is inversely proportiona l to r (e.g., iPotekhinl 
1999; lYakovlev fc Kaminkerlll99l . 

4.3. Total photon luminosities 

In order to calculate the cooling curves, one needs the 
total photon luminosity L 1 as a function of Tb. In the 
magnetic case, T 7 is given by the average over the sur- 
face, Eq. ©. 

Figure |H1 displays the photon luminosity versus B for 
two selected values of T b (3 x 10 7 and 3 x 10 s K) and 
four selected values of AM. The dependence of L 7 on 
B is complicated. At T b = 3 x 10 8 K (in a warm NS) 
and not very strong fields, the equatorial reduction of the 
heat transport clearly dominates, and the NS luminosity 
is lower than at B = 0. For higher B, the polar enhance- 
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Fig. 7. — T s vs. Tt, for the canonical NS with fully accreted (solid line and dots), iron (dot-dashed lines and triangles) and partly accreted 
envelope at AM/M = 10 -12 (dashed lines and empty squares), for 9 = (upper lines and symbols) and 9 = 90° (lower ones). The lines 
show the fit, and the symbols show the numerical results. Left, middle, and right panels: B = 10 11 G, B = 10 13 G, and B = 10 15 G, 
respectively. 
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Fig. 8. — Photon surface luminosity (redshifted as detected by 
a distant observer, left vertical axis) or redshifted effective sur- 
face temperature (right vertical axis) of a canonical NS with a 
dipole magnetic field, for two values of T\, and four models of 
the heat- blanketing envelope (accreted mass AM = 0, 10 — 12 Mq, 
10 -10 Mq, or 10 — 7 Mq) versus magnetic field strength at the mag- 
netic pole. 



ment of the heat transport becomes more important, and 
the magnetic field increases the photon luminosity. At 
lb = 3 x 10 7 K (in a much colder NS) the effect of mag- 
netic fields B < 10 13 G is very weak while the increase 
of L 7 by stronger fields is much more pronounced. This 
is because the electron contribution into conduction be- 
comes lower (and the radiative contribution higher) in a 
cold plasma. Accordingly, the equatorial decrease of the 
heat transport (associated with electron conduction) is 
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Fig. 9. — Local surface temperature (in 10 6 K) vs. cosine of 
the polar angle \ f° r a NS model with M = 1.4 Mq, R = 10 km, 
T b = 2 X 10 s K, B = (horizontal lines) and the dipole field with 
Bp = 10 G (curved lines), for nonaccreted (dashed lines) and 
accreted (solid lines) envelopes. 



weaker. On the other hand, the quantum effects increase 
both, the electron and radiative thermal conductivities, 
and are more pronounced in a colder plasma. Generally, 
the photon luminosity is not a simple function of B, be- 
cause it is affected by a number of factors. Radiative 
opacities, longitudinal electron conductivities, position 
of the radiative surface all depend on B and affect T s . 
In Fig. |S| the variation of L 1 for B p < 10 14 G does not 
exceed a factor of 2.5. 

A higher B causes a more significant increase of the 
photon luminosity. However, in this case the effect of an 
accreted envelope is weaker than at B = 0. The origin 
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Fig. 10. — Cooling ol nonsuperfluid NSs with M = 1.3 Mq (solid lines, slow neutrino emission) and 1.5 Mq (dashed lines, fast neutrino 
emission) assuming EOS B in the NS cores. Left panel: nonmagnetic (NM) NSs with different amount AM of light elements in the heat 
blanketing envelopes (the values of log AM/M are shown near the curves); thick curves refer to nonaccreted (Fe) envelopes. Middle panel: 
nonaccreted envelopes with dipole surface magnetic fields (the curves are marked with logB, where B = B p is the field at the pole); thick 
lines refer to B = 0. Right panel: same as in the middle panel but with a fully accreted envelope. 



of this weakening is explained in Fig. |3 which illustrates 
the dependence of the local effective temperature on the 
magnetic latitude, specified by the polar angle x? Eq. 
HU, for a canonical NS with T b = 2 x 10 8 K. The dashed 
and solid lines refer to the nonaccreted and fully accreted 
envelopes, respectively. The horizontal lines refer to the 
case of B = in which the accreted envelope increases 
T e g by a factor of 1.73. The curved lines show the case 
of B p = 10 15 G. Since the effects of the magnetic field 
weaken with increasing temperature (Paper II), the in- 
crease of T s at the magnetic pole for the hotter accreted 
envelope is smaller than the analogous increase for the 
cooler iron envelope. Hence the equatorial band, where 
T s decreases, is broader for the accreted envelope, lead- 
ing to an additional compensation of the total luminosity 
increase. 

5. COOLING 

Let us outline the effects of accreted envelopes and 
surface magnetic fields on NS cooling. We use the 
sa me nonisothermal, g eneral relativistic cooling code as 
in iGnedin et alJ 1)200 1|) . but incorporated the above ef- 
fects of accreted envelopes and surface magnetic fields. 
In particular, we have shifted the inner boundary of the 
heat blanketing envelope to the neutron-drip density (see 

a. 

For simplicity, we assume that the NS cores are com- 
posed of neutrons, protons, and electrons. We adopt two 
model EOSs of this ma tter, EOS A and EOS B, based on 
the E OSs proposed bv lPrakash. Ainsworth. fe Lattimerl 
1 1988) and described, for instan ce, in lYakovl ev et al.l 
(2002b). EOS A is model I of iPrakash et alJ <|1988ft 
with the compression modulus of saturated nuclear mat- 
ter K = 240 MeV. EOS B corresponds to K = 180 
MeV and to the simplified form of the symmetry energy 
proposed by iPaee fe Applegatd l)1992|) . The NS models 
based on EOSs A a nd B are described, for instance, by 
IGnedin et alJ l)2001l) . EOS A is somewhat stiffer, and 
yields the maximum NS mass 1.98 Mq, while EOS B 
yields the maximum mass of 1.73 Mq. Both EOSs al- 



low the powerful direct Urea process of neutrino emis- 
sion to operate at sufficiently high densities p > (with 
p D = 7.85 x 10 14 g cm" 3 and 1.298 x 10 15 g cm' 3 for 
EOSs A and B, respectively). We consider the NS models 
of two masses, 1.3 and 1.5 Mq. The parameters of these 
models are listed in Tabled Pc is the NS central density, 
M crust the crust mass, Ai? crust the crust thickness (de- 
fined as Ai? crust = R— i? cc , -Rcc being the circumferential 
radius of the crust-core interface) , Md is the mass of the 
inner core (if available) where direct Urea process oper- 
ates, and i?D is the circumferential radius of this core. 
The central densities of the low-mass NSs, 1.3 M©, are 
smaller than p-Q for both EOSs. They give an example of 
slow cooling. The central densities of 1.5 Mq NSs exceed 
Pd, i.e., these models give us an example of fast cooling. 

5.1. Overall effects 

Figure shows the effects of the surface magnetic 
fields and accreted envelopes on the cooling of 1.3 (solid 
lines) and 1.5 Mq (dashed lines) nonsuperfluid NS mod- 
els with EOS B in the core. 

The left panel of Fig. ^| illustrates the effects of ac- 
creted envelopes in nonmagnetized NSs. We present the 
cooling curves for some values of AM, the mass of rel- 
atively light elements (H, He, C, and O) in the heat- 
blanketing envelopes. The curves for nonaccreted (Fe) 
envelopes are plotted by thick lines. The fraction of 
accreted mass AM/M varies from (nonaccreted en- 
velopes) to ~ 10~ 7 (fully accreted envelopes); a further 
increase of AM is li mited by the pycnon uclear burning 
of light elements (cf. iPotekhin et alJll997[) . 

During the first 50 years after the birth, 1.3 Mq and 
1.5 Mq NSs have nearly the same surface temperatures 
since the surface is thermally decoup led from the stel- 
lar interiors (e.g.. IGnedin et aTll200lD . Later, after the 
thermalization, the direct Urea process in the 1.5 Mq NS 
makes this NS much colder. The change of slopes of the 
cooling curves at t ~ 10 5 yr reflects transition from the 
neutrino to the photon cooling stage. At the neutrino 
stage, the internal stellar temperature is ruled by the 
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Table 1. Neutron star models 



EOS 


M 
(Mq) 


R 

(km) 


Pc 

(10 14 g cm- 3 ) 


-^crust 

(Mq) 


A^crust 

(km) 


Md 

(Mq) 


Ru 
(km) 


A 


1.3 
1.5 


13.04 
12.81 


7.44 
9.00 


0.057 
0.049 


1.58 
1.26 


0.137 


4.27 


B 


1.3 
1.5 


11.86 
11.38 


10.70 
14.20 


0.039 
0.028 


1.26 
0.93 


0.065 


2.84 



neutrino emission and is thus independent of the ther- 
mal insulation of the blanketing envelope. The surface 
photon emission is determined by the T^-T e f{ relation. 
Since the accreted envelopes of not too cold NSs are more 
heat transparent, the surface temperature of an accreted 
star is noticeably higher than that of a nonaccreted one. 
One can see that even a very small fraction of accreted 
matter, such as AM/M ~ 10~ 13 , can change apprecia- 
bly the thermal history of the star. The colder the star, 
the smaller the fraction of accreted material which yields 
the same cooling curve as the fully accreted blanketing 
envelope. This effect is more pronounced for the fast 
cooling. 

At t > 10 5 yr, a star enters the photon cooling stage. 
Since the accreted blanketing envelopes have lower ther- 
mal insulation, the NSs with such envelopes cool faster 
at the photon stage than the nonaccreted NSs (Fig. HU|) . 
Thus, the light elements make the opposite effects on the 
surface temperature at the neutrino and photon cooling 
stages. This reversal of the effect while passing from one 
stage to the other is very well known and quite natu- 
ral. Similar results have been obtained in Paper I, but 
our new T e g — Tt> relation slightly weakens the effect of 
accreted envelopes . 

The middle panel of Fig.lTUldisDlavs the effect of dipole 
magnetic field on the cooling of NSs with nonaccreted en- 
velopes. We present the cooling curves for several mag- 
netic field strengths at the magnetic poles (numbers next 
to the curves) up to B p = 10 16 G. The cooling curves of 
nonmagnetic NSs are plotted by thick lines. For simplic- 
ity, the magnetic field is treated as fixed (nonevolved) . 
The thermal state of the stellar interior is almost inde- 
pendent of the magnetic field in the NS envelope at the 
neutrino cooling stage, but is affected by the magnetic 
field later, at the photon cooling stage. On the contrary, 
the surface temperature is always affected by the mag- 
netic field. 

The dipole field B p < 10 13 G makes the blanket- 
ing envelope of a warm (1.3 Mq) NS overall less heat- 
transparent ( H4.3p . This lowers T e fi at the neutrino cool- 
ing stage and slows down cooling at the photon cooling 
stage. The dipole field B p 3> 10 13 G makes the blanket- 
ing envelope overall more heat transparent, increasing 
T e g at the neutrino cooling stage and accelerating the 
cooling at the photon stage. The field B p ~ 10 13 G has 
almost no effect on the NS cooling. The presented results 
are in satisfactory agreement with our previous studies 
(Paper II). 

The rapid cooling of colder magnetized (1.5 Mq) NSs 
is somewhat different. A strong magnetic field makes 
the heat-blanketin g en velopes of these NSs overall more 
heat transparent ( M4.3[) . i.e., increases T c g at the neutrino 



cooling stage and decreases it at the photon cooling stage. 
The fields B p < 10 13 G have almost no effect on T g. On 
the contrary, the effects of higher fields, B p > 10 13 G, 
are much stronger than for slowly cooling NSs. 

The right panel of Fig. EH presents the cooling curves 
of NSs with fully accreted envelopes and the same dipole 
magnetic fields as in the middle panel. For a NS with 
B p < 10 15 G at the neutrino cooling stage the effect 
of the accreted envelope is stronger than the effect of 
the magnetic field. For higher B p , the magnetic effect 
dominates; the accreted envelope produces a rather weak 
additional rise of T c g . 

The main outcome of these studies is that even ultra- 
high magnetic fields cannot change the average surface 
temperatures of young and warm NSs as appreciably as 
an accreted envelope can, although the distribution of 
the local surface temperature over the surface of a mag- 
netized NS can be strongly nonuniform. 

5.2. Very slowly cooling low-mass NSs 

As demonstrated rec e ntly by 

I Yakovlev. Kaminker. & Gnedinl ll2001aD : I Yakovlev et alJ 
l|2002b|) and lKaminker et alJ l)2002j) . the effects of strong 
magnetic fields and accreted envelopes are especially 
important in low-mass NSs (where the direct Urea 
process is forbidden) with a strong proton superfluidity 
in their cores. Such a superfluidity (with typical values 
of proton critical temperature T cp (p) >5x 10 9 K) fully 
suppresses the modified Urea process which otherwise 
would be the main neutrino emission mechanism in 
low-mass NSs. The neutrino emission is then generated 
in the reactions of neutron-neutron scattering; it is 
not suppressed by a proton superfluidity and becomes 
dominant. 

The low-mass NSs with strongly superfluid protons 
have thus very low neutrino luminosity and form a special 
class of very slow-cooling stars. Their models are useful 
to interpret the observat ions of isolated NSs hottest for 
thei r ages. According to [ Yakovlev et alJ l)2001al [2002a) 
and iKaminker et alJ l|2002fl . there are two NSs of such 
a type among several isolated NSs whose thermal radi- 
ation has been detected. One of them, RX J0822— 43, 
is young, while the other, PSR B1055— 52, is much 
older. The third possible candid ate, RX J1856-3754, 
was moved fYakovlev et alJl2002aj) from the class of very 
slow-cooling NSs to the class of faster coo lers after the 
revision of its age {Walter fc Lattimerl2002fl . In this sub- 
section, we focus on the interpretation of the observa- 
tions of RX J0822-43 and PSR B1055-52, taking into 
account the effects of magnetic fields and accreted en- 
velopes. The latter effects are less important for the in- 
terpretation of observations of colder isolated NSs (see, 
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Fig. 11. — Energy gaps (left vertical axis) or critical temper- 
ature (right vertical axis) for various models of crustal neutron 
superfluidities (see the text) versus neutron Fermi wave number; 
the vertical dotted line marks the crust-core interface. 



e.g. , iKaminker et al.ll2002l and lYakovlev et~aTll2002bl for 

a recent comparison of observations and cooling theories 
of NSs with nucleon cores). 

Our analysis is illustrated in Figs. Illlandll2l The ob- 
servational limits T$ = (1.6-1.9) MK (at the 90% con- 
fidence level) for RX J0822-43 (F ig. IT2jl are taken from 
IZavlin. Trumper. fe Pavlovl l|1999|) . They are obtained 
by fitting the observed spectrum with the hydrogen at- 
mosphere model; such a model is more appropriate for 
this source than the blackbody model of thermal emis- 
sion. In contrast, it seems more reasonable to fit the 
spectrum of PSR B1055— 52 with a model containing a 
blackbody thermal component f e.g. . [Pavlov et alH 2002). 
The results of different groups are not fully consistent 
because of the complexity of observations and their inter- 
pretation. For instance, according to G. G. Pavlov (2002, 
private communication) the best-fit value obtained from 
the recent Chandra observations of PSR B1055— 52 is 
« 60 eV. Attributing a typical error bar (at the 90% 
confidence level) derived from the observations of the 
same source somewhat earlier (G. G. Pavlov & M. Teter 
2002, private communication) we have T£g = (60 ± 6) 
eV. An analysis of the recent BeppoSAX ob servations of 
PSR B1055-52 gives T e g = (75 ± 6) eV l|Mineo et all 
2002) at the 68% confidence level. We adopt, somewhat 
arbitrarily, a wide interval from 54 eV (the minimum 
value from the Chandra data) to 81 eV (the maximum 
value from the BeppoSAX data). We expect that this 
interval reflects the actual uncertainty of our knowledge 
of for PSR B1055-52. 

The NS ages are also known with some uncertainty. 
For RX J0822-43, we take the age range t = 2 
5 kyr (as can be deduced, e.g., from a discussion in 
lArendt. Dwek fc Petrel 11988ft centered at t = 3.7 kyr 
fWinkler et el .111 9881) . For PSR B1055-52, we adopt the 
standard spin-down age of 530 kyr and assume that it is 
uncertain within a factor of 2. 

For our analysis, we take a 1.3 Mq NS model with 
EOS A in the core. We use model lp of sufficiently 
strong proton superfluidity in the core (in notations of 



IKaminker et al.ll2002ft and assume a weak (triplet-state) 
neutron pairing in the core with the maximum critical 
temperature lower than 2 x 10 s K. This weak neutron 
superfluidity has no effect on NS cooling and can be 
neglected in the cooling simulations. A stronger neu- 
tron superfluidity in the NS core would initiate a pow- 
erful neutrino emission due to Cooper pairing of neu- 
trons accelerating NS cooling in disagreem ent with the 
ob servations (e.g.. lYakovlev et alJ l2002albft . As shown 
by lYakovlev et al.l l)2001aft . the cooling of low-mass NSs 
is rather insensitive to the model EOS in the NS core, to 
the NS mass (as long as the mass is sufficiently low to 
avoid fast neutrino cooling in the NS core), and to the 
model of proton superfluidity in the NS core (as long as 
the proton critical temperature is > (2 — 3) x 10 9 K over 
the core to suppress the modified Urea process). There- 
fore, model lp of proton superfluid is used just as an ex- 
ample, and the cooling curves are actually independent 
of the features of strong proton superfluid. Let us add 
that it is likely that the cores of low-mass NSs consist of 
nucleons and electrons (with forbidden direct Urea pro- 
cess) . Muons may also be present there, bu t have almost 
no ef fect on the cooling of low-mass NSs l)Beiger et alJ 
l20?)l . Thus, our cooling scenarios are not related to the 
models of matter of essentially supranuclear density in 
the inner cores of massive NSs where the composition 
may be exotic (e.g., includes pion or kaon condensates, 
or quark matter). 

We see that the models of low-mass cooling NSs are 
sufficiently robust against uncertainties in the physics 
of matter in the NS cores. Accordingly, the cooling of 
this special class of NSs (contrary to the cooling of other 
NS models) is especially sensitive to the properties of the 
NS crust. It is mainly regulated by the effects of (i) 
accreted matter and (ii) surface magnetic fields in the 
heat-blanketing envelopes, as well as by the effects of (Hi) 
singlet-state neutron superfluidity in the inner NS crusts. 
All these effects are of comparable strength. They are an- 
alyzed below and illustrated in Figs. II II and 1121 with the 
aim to interpret the observations of RX J0822— 43 and 
PSR B1055-52. 

Two dot-and-dashed curves marked as noSF and pSF 
in each panel of Fig. 1121 show the cooling of a nonsuper- 
fluid NS and a NS with strong proton superfluidity in 
the core (the effects of magnetic fields and accreted en- 
velopes are neglected). We see that the proton super- 
fluidity, indeed, delays the cooling (by suppressing the 
modified Urea process) and makes the cooling curves con- 
sistent wit h the observations of RX J0822— 43 an d PSR 
B1055-52 (|Kaminker. Haensel. fc Yakovlevll200l . 

However, our "successful" pSF cooling curve neglects 
the presence of neutron superfluidity in the NS crust. 
This superfluidity is predicted by microscopic theories 
(e.g.. lLombardo fc Schufzel2001l) although the superfluid 
gaps are very model dependent. The superfluidity initi- 
ates neutrino emission due to Cooper pairing of crustal 
neutrons whic h may noticeably accele rate the cooling of 
low-mass NSs ( Yakovlc v et alJfeOOlaj) . In order to illus- 
trate this effect, we consider several models of the crustal 
superfluidity. Figure ITT1 shows the dependence of the su- 
perfluid gap A (left vertical axis) or associated critical 
temperature T cn (in units of 10 9 K, right vertical axis) on 
neutron Fermi wave numb er fcFn (as a measure of den - 
sity) for six models (from iLombardo fc Schulzel 12001). 
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Fig. 12. — Evolution of redshifted effective surface temperature (T£g) or photon luminosity (£2°) of a low-mass (1.3 Mq) NS with 
EOS A confronted with observations of RX J0822— 43 and PSR B1055— 52. Dot-and-dashed curves: nonmagnetized nonaccreted NS model 
without any superfluidity (noSF) or with a strong proton core supcrfiuid alon e (pSF). Other curves are obtained including the effects of 
proton core superfluidity and a model of crustal neutron pairing from Fig. 1111 fBCS or S: left panel; W or A: middle panel; C86 or C93: 
right panel). Models Sa7, Wa8, and C93a9 include also the effects of accreted envelopes with AM/M = 10 — 7 , 10 — 8 , and 10 -9 , respectively. 
Models Sb, Wb, and C93b refer to nonaccreted envelopes with dipole magnetic field B p = 10 15 G. 



Vertical dotted line indicates approximate position of the 
crust-core interface (assumed to be at 1.5 x 10 14 g cm~ 3 
in our cooling models). 

Model BCS is the basic model of singlet-state 
neutron pairing calculated under the simplified as- 
sumption of purely in-vacuum neutron-neutron in- 
teraction (neglecting mediu m polarization e ffects) . 
Other five mode ls - C86 llChen et alJ llQSl. C93 
llChen et all 11993ft. A llAinsworth. Wambach fc Pines! 
I- 1 9891). W iWarnhach. Ainswortn. fc Pined 11 99.3ft. and S 
(Schulze et al. 199Q) - include medium polarization ef- 
fects which weaken the strength of neutron pairing. 

In Fig. ^1 in addition to the noSF and pSF curves, 
we show six cooling curves (BCS, S, A, W, C86, and 
C93) calculated adopting proton superfluidity in the NS 
core and one of the models of neutron superfluidity in the 
crust from Fig. II II (neglecting surface magnetic fields and 
accreted envelopes) . One can see that the crustal super- 
fluidity, indeed, accelerates the cooling and complicates 
the interpretation of the observations of RX J0822— 43 
and PSR B1055— 52. The new six curves are naturally 
divided into three pairs shown in three panels of Fig. 1121 
(left panel: BCS and S; middle panel: A and W; right 
panel C86 and C93). The curves within each pair are 
very close, while the pairs differ from one another. As 
clear from Fig. ^] the neutron superfluid gaps for any 
pair, although different, have one common property: the 
same maximum density of superfluidity disappearance. 
This maximum density limits the density of the outer NS 
layer where the neutrino emission due to the singlet-state 
pairing of neutrons operates and accelerates the cooling. 
For instance, neutron superfluids BCS and S extend to 
higher densities (penetrate into the NS core); Cooper 
pairing neutrino emission is generated from an extended 
layer and produces the most dramatic effect: the cool- 
ing curves go much lower than the pSF curve, strongly 
violating the interpretation of RX J0822— 43 and almost 



violating the interpretation of PSR B1055— 52. Superflu- 
ids A and W die at lower p (at the crust-core interface), 
Cooper-pairing neutrinos are emitted from a smaller vol- 
ume, and the cooling curves go higher, closer to the pSF 
curve. They do not explain RX J0822— 43 but seem to 
explain PSR B1055-52. Finally, superfluids C86 and 
C93 disappear even at lower p, far before the crust-core 
interface; the emission volume gets smaller, and the cool- 
ing curves shift even closer to the pSF curve simplifying 
the interpretation of RX J0822— 43 and easily explaining 
PSR B1055-52. 

All NSs should have the same EOS and superfluids in 
their interiors but may have different magnetic fields and 
accreted envelopes. The effects of latter factors are also 
illustrated in Fig.^J As discussed above, the cooling his- 
tories of NSs with the superfluid model BCS are almost 
the same as with S, with model A are the same as with 
W, and with model C86 are the same as with C93. There- 
fore, we do not consider the effect of magnetic fields and 
accreted envelopes on the models C86, A and BCS. Let 
us remind, that magnetic fields and accreted envelopes 
have opposite effects on the thermal states of NSs at the 
neutrino and photon cooling stages. PSR B1055— 52 is 
just passing from one cooling stage to the other and has 
no superstrong magnetic field. It is not expected to pos- 
sess an extended accreted envelope. Thus, the effects of 
magnetic fields and accreted envelopes on the evolution 
of this pulsar are thought to be minor. 

The superfluid model S (or BCS) without any magnetic 
field or accreted envelope is only marginally consistent 
with the observations of PSR B1055— 52 (left panel of 
Fig. I12|l . If, however, we accept this model, then we can 
explain RX J0822— 43 by switching on the effects of the 
accreted envelopes or surface magnetic fields. Curve Sb 
on the left panel of Fig. El is calculated adopting proton 
superfluid in the NS core, superfluid S in the crust, and 
the dipole magnetic field with B p = 10 15 G (the mag- 
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netar hypothesis). Curve Sa7 is obtained for the same 
crustal superfiuid, but for B = and for AM/M = 1(T 7 
of accreted material on the NS surface. We see that 
both curves, Sb and Sa7, are consistent with the ob- 
servations of RX J0822— 43. The magnetic field slightly 
above B p = 10 15 G would further improve the agreement 
of the theory and observations. 

The next superfiuid model W (or A) is acceptable to 
interpret the observations of PSR B1055-52 (Fig. El 
middle panel) . We can adopt it and switch on the effects 
of the magnetic field or accreted envelope to interpret RX 
J0822— 43. Curve Wb on the right panel corresponds to 
the crustal superfiuid W and the dipole magnetic field 
with Bp — 10 15 G. Curve Wa8 is calculated for the same 
crustal superfiuid, B = 0, but for AM/M = 1CT 8 . Both 
curves, Wb and Wa8, are seen to be consistent with the 
observations of RX J0822-43. 

Finally, model C93 (or C86) of crustal superfiuid is in 
reasonable agreement with the data on PSR B1055— 52 
(Fig- El right panel). Adding the effect of the surface 
magnetic field (B p = 10 15 G, curve C93b) or the accreted 
envelope (AM/M = 10~ 9 , curve C93a9) we can easily 
explain the observations of RX J0822— 43. 

Thus all the models of crustal superfluidity are cur- 
rently consistent with the observations (although models 
BCS and S seem to be less likely for explaining PSR 
B1055— 52). This is a consequence of wide observational 
error bars of for PSR B1055-52 (see above). We 
expect to constrain the models of crustal superfluidity 
after better determination of in the future observa- 
tions of PSR B1055— 52. Afterwards it will be possible 
to constrain the surface magnetic fields and the mass 
of the accreted envelope from the observations of RX 
J0822-43. 

It is clear from Figs. El an d El that old and warm 
NSs like PSR B1055— 52 cannot possess bulky accretion 
envelopes which would operate as efficient coolers. This 
conclusion is in line with the observations of thermal ra- 
diatio n from old sources : their spectra are better fitted 
(e.g., iPavlov et al.l 12002^1 with the blackbody model of 
thermal emission (suitable for nonaccreted matter) than 
with the hydrogen atmosphere models. This conclusion is 
also in line with theoretical studies of lChang fc Bildstenl 
( 2003) who show that light elements can be burnt in old 
(t > 10 5 yr) and warm NSs by diffusive nuclear burning. 
In this connection, it would be interesting to construct 
the models of cooling NSs incorporating the effects of 
diffusive burning and associated thinning of the light- 
element envelopes in time. Cooling of magnetars is ex- 
pected to be accompanied by magnetic field decay which 
should also be taken into account in advanced cooling 
simulations. 

6. CONCLUSIONS 

We have studied the thermal structure of a heat- 
blanketing envelope of a NS with a strong magnetic field 



and arbitrary amount of light-element (accreted) ma- 
terial. We have calculated and fitted by an analytic 
expression the relation between the neutron-star inter- 
nal and surface temperatures in a local element of the 
heat-blanketing envelope as a function of magnetic field 
strength and geometry, the mass of accreted material, 
and the surface gravity. We have performed numerical 
simulations of cooling of NSs with dipole magnetic fields 
and accreted envelopes. We have considered slow and 
fast cooling regimes but mainly focused on a very slow 
cooling of low-mass NSs with strong proton superfluid- 
ity in their cores. These NSs form a special class of NSs 
whose thermal history is rather insensitive to the physics 
of matter in the stellar cores but is mainly regulated by 
the magnetic field strength and the amount of accreted 
material in the heat-blanketing envelopes, as well as by 
the singlet-state pairing of neutrons in the inner stel- 
lar crusts. We show that these cooling regulators are 
important for explaining the observations of thermal ra- 
diation from isolated NSs warmest for their ages, RX 
J0822-43 and PSR B1055-52. Our analysis indicates 
that all realistic microscopic models of the singlet-state 
neutron pairing are currently consistent with the obser- 
vations of PSR B1055— 52. Adopting these models and 
tuning the strength of the magnetic field and/or the mass 
of accreted material we can explain also the observations 
of RX J0822— 43. We expect that such an interpreta- 
tion will be refined after new, high-quality observations 
of RX J0822-43 and PSR B1055-52 will appear or new 
very slowly cooling NSs will be discovered. In particular, 
this will allow one to discriminate between the models of 
crustal superfluidity and determine the depth of super- 
fluidity disappearance in high-density matter. It will also 
be possible to constrain the surface magnetic fields and 
the mass of accreted envelopes of these objects. Our re- 
sults will also be useful for constructing advanced models 
of cooling NSs taking into account the evolution of strong 
surface magnetic fields and the mass of light-element en- 
velopes (e.g., under the action of diffusive nuclear burn- 
ing). 

We believe that the T^-T s relation obtained in this 
paper is reliable at B < 10 14 G. The results presented 
for higher fields are rather indicative, because the dense 
plasma effects on the heat conduction by photons near 
the bottom of NS photosphere at the superstrong fields 
have not yet been explored in detail. 
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APPENDIX 

FITTING FORMULAS FOR THE T h -T s RELATION 

Let T b9 = T b /10 9 K, T s6 = T s /10 6 K, and r) = g\ A AM /M, where 314 = .9/IO 14 cm s" 2 and AM is the mass of 
accreted light elements from H to O. We assume that ij < 10 -6 ; at larger r\ the light elements would undergo efficient 
pycnonuclear burning. The envelope is fully accreted if ry = 10~ 6 , and nonaccreted if 77 = 0. Let T a § and Tp e 6 denote 
T S 6 for the fully accreted and nonaccreted envelopes, respectively. 
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Table 2. Parameters of Eqs. IA5I . IA6I 



envelope 


n 


1 




2 


3 


4 


5 


6 


7 


8 


9 


10 


11 




iron 


On 
b„ 


1.76E- 
159 


4 


0.038 
270 


1.5 
172 


0.0132 
110 


0.620 
0.363 


0.318 
0.181 


2.3E-9 
0.50 


3 

0.619 


0.160 


21 


4.7E- 


-5 


accreted 


dn 


4.50E- 
172 


3 


0.055 
155 


2.0 

383 


0.0595 
94 


0.328 
0.383 


0.237 
0.367 


6.8E-7 
2.28 


2 

1.690 


0.113 


163 


3.4En 


-5 



First consider the case of B = 0. In this case, the fit given by Eqs. (A6) and (A7) of Paper I for the nonaccreted 
envelope remains valid: 



TL = 9u [(7C) 2 - 25 
For the fully accreted envelope, we have 



(0.33C) 125 ], C=7b9-10- 3 ff 1 Tv / 77^ 



rri4 



g u (18-1 T h9 ) 2A2 {0.447 + 0.075 (log 10 T b )/[1 + (6.2 T b9 ) 4 ]} + 3.2 T^ 67 Tp c6 /(l + 3.2 T b g 



.67\ 



(Al) 
(A2) 



This fit differs from the equations in Appendix C of Paper I by the correction factor in curly brackets. The difference 
is most appreciable at Tb > 10 8 K, as explained in H3.ll (cl. Fig.EJ). 
Next consider the fully accreted and nonaccreted magnetized envelopes. In both cases, we can write 



TFeG,a6{B) = TFo6,a6(0)<%', 

where the magnetic correction factor X depends on B, Tb, and 9. It can be fitted by 

1/a 



X - 



X\\ 



(*ii 



cos 2 6 + X? sin 2 6 



4+^XjXj { forT Fc6 , 
(2 + X ± /X l{ ) 2 forT a6 , 



13 12 



T^ + a^l (l+a 7 B 12 /T$) 



3.7 + (a w + a n B 12 3/2 ) T x 



^2 

l b9, 



1 + h B 12 /(l + b 2 T°l) 



1/2 



l + b 3 B 12 /(l + b 4 T**) 



/?= (l + 6 B 7*S 



(A3) 

(A4) 
(A5) 

(A6) 



with the parameters a, and bi given in Table 
Finally, the surface temperature of a partly accreted envelope is approximately reproduced by the interpolation: 



r s = [7^ 6 + (i-7)^ c6 ] 1/4 , 

7 = [l + 3.8 (0.10 9 ] _1 [l + 0.171 C 7/2 T\ 



l:>9 



e = -log 10 (10 6 77). 



(A7) 
(A8) 



These fitting formulas have been checked against calculations for input parameters restricted by the conditions 
6.5 < log 10 Tb/K < 9.5, log 10 T s /K > 5.3, and 10 < log 10 B/G < 16. The numerical values of T s are reproduced with 
root-mean-square residuals of (3-5)% (< 2% for fully accreted envelopes at — 0) and maximum deviations within 
w 20% (within 13% for 8 = 0, within 12% for fully accreted or nonaccreted envelopes at any 9, and within 5.2% for 
fully accreted envelopes at 9 = 0). 

We emphasize that our numerical results, and hence the fitting formulas given here, are uncertain at superstrong 
fields (B > 10 14 G) because of the dense plasma effects discussed in M3.3I fsee the right panel of Fig.Q. A thorough 
study of the radiative heat conduction at lu < u> p \ is needed to obtain a reliable T^-Tg relation at such field strengths. 
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